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(Dated: February 27, 2006) 

We present the measurement of 7i = B{t — > Wb)/B{t —^ Wq) in pp collisions at ^/s = 1.96 TeV, 
using 230 pb^^ of data collected by the D0 experiment at the Fermilab Tevatron Collider. We fit 
simultaneously R and the number (iN/jf) of selected top quark pairs (tt), to the number of identified 
6-quark jets in events with one electron or one muon, three or more jets, and high transverse energy 
imbalance. To improve sensitivity, kincmatical properties of events with no identified 6-quark jets 
are included in the fit. We measure R = I.OSIIq'I^ (stat + syst), in good agreement with the standard 
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model. We set lower limits of i? > 0.61 and \Vtb\ > 0.78 at 95% confidence level. 



PACS numbers: 12.15.Hh, 14.65.Ha 

Within the standard model (SM), the top quark de- 
cays 99.8% of the time to a.W boson and a b quark, with 
the ratio i?= B{t Wb)/B{t Wq) (here q refers to 
d, s, or b quarks) expressible in terms of the Cabbibo- 
Kobayashi-Maskawa (CKM) matrix elements |1] i? = 

|Vt^P-i- | V't''h+|Vtj|^ • '^^^ unitarity of the CKM matrix and 
experimental constraints on its elements yield the SM 
prediction 0.9980 < R < 0.9984 at the 90% C.L. Nev- 
ertheless, a fourth generation of quarks or non-SM pro- 
cesses in the production or decay of the top quark could 
lead to significant deviations from the SM. So far, mea- 
surements of R by the CDF collaboration ^ have not 
established a deviation of R from unity. 

In the present analysis, we assume that the top quark 
decays into a W boson, but that the associated quark 
can be d, s, or b. Lepton -I- jets final states arise in tt 
when one W boson decays leptonically and the other into 
a qq' pair. About 6% of the signal arises from tt events in 
which both W bosons decay leptonically, but one charged 
lepton is not reconstructed, while additional jets are pro- 
duced by initial or final state radiation. In this Letter, 
we report the measurement of R in the lepton (electron 
or muon) + jets channel {£ -\- jets). The lepton can come 
either from a direct W decay or from W ^ t ^ e//Lt. We 
use 6-jet identification (6-tagging) techniques, exploiting 
the long lifetime of B hadrons, to separate ti events from 
the background processes. The data were collected by the 
D0 experiment from August 2002 through March 2004, 
and correspond to an integrated luminosity of 230 pb~^. 

The D0 detector incorporates a tracking system, 
calorimeters, and a muon spectrometer j^]. The tracking 
system is made up of a silicon micro-strip tracker (SMT) 
and a central fiber tracker (CFT), located inside a 2 T 
superconducting solenoid. The tracking system provides 
efficient chargedparticle detection in the pseudorapidity 
region \r]\ < 3 The SMT strip pitch of 50-80 ^m 
allows a precise determination of the primary interaction 
vertex (PV) and an accurate measurement of the impact 
parameter of a track relative to the PV • These are key 
components of the lifetime-based ^-tagging algorithms. 
The PV is required to be within the fiducial region of the 
SMT and to contain at least three tracks. The calorime- 
ter consists of a barrel section covering |ry| < 1.1, and 
two end-caps extending the coverage to \r]\ « 4.2. The 
muon spectrometer surrounds the calorimeter and con- 
sists of three layers of drift chambers and several layers 
of scintillators |8j. A 1.8 T iron toroidal magnet is lo- 
cated outside the innermost layer of the muon system. 
The luminosity is calculated from the rate of pp inelastic 
collisions, detected by two arrays of scintillation counters 
mounted close to the beam-pipe on the front surfaces of 



the calorimeter end-caps. 

We select data in the electron and muon decay chan- 
nels by requiring an isolated electron with px > 20 GeV 
and \t]\ < 1.1, or an isolated muon with pT > 20 GeV 
and |ry| < 2.0. The lepton isolation criteria are based on 
calorimeter and tracking information. More details on 
lepton identification and trigger requirements are avail- 
able in Ref. 0. In both channels, we require the miss- 
ing transverse energy {^t) to exceed 20 GeV and not 
be colinear with the direction of the lepton projected 
on the transverse plane. The candidate events must be 
accompanied by jets with pT > 15 GeV and rapidity 
\y\ < 2.5 0. Jets are defined using a cone algorithm 
with radius A7^ = 0.5 [l3|. 

We use a secondary vertex tagging (SVT) algorithm to 
reconstruct displaced vertices produced by the decay of B 
hadrons inside jets. Secondary vertices are reconstructed 
from two or more tracks satisfying: px > 1 GeV, > 1 hits 
in the SMT detector, and impact parameter significance 
dca/Sd,a > 3.5 Tracks identified as arising from Ka 
or A decays or from 7 conversions are not used. If the 
secondary vertex reconstructed within a jet has a decay- 
length significance L^y/Sh^y > 7 |0, the jet is defined 
as 5-tagged. Events with exactly 1 (> 2) 6-tagged jets 
are referred to as 1-tag (2-tag) events. Events with no 
6-tagged jets are referred to as 0-tag events. A prediction 
for the number of background events and the fractions of 
tt events in the 0, 1, and 2-tag samples require the prob- 
abilities for different types of jets (5-, c-, and light-quark 
jets) to be 5-tagged. The calculation of these probabil- 
ities is presented in Ref. We fit simultaneously R 
and the total number of ti events in the 0, 1, and 2- 
tag samples (iVtt) to the number of observed 1-tag and 
2-tag events, and, in 0-tag events, to the shape of a dis- 
criminant variable V that exploits kinematic differences 
between the backgrounds and the ti signal. 

The main background in this analysis is from the pro- 
duction of leptonically decaying W bosons produced in 
association with jets (Vl^-l-jets). Most of the jets ac- 
companying the W boson originate from u, d, and s 
quarks and gluons (VF -I- light jets). Between 2% and 14% 
of M^-|-jets events contain heavy-ffavor jets, arising from 
gluon splitting into bb or cc {Wbb or Wcc, respectively). 
About 5% of the M^-|-jets events contain a single c quark 
that originates from VF-boson radiation from an s quark 
in the proton or anti-proton sea (s —^ Wc). A sizable 
background arises from strong production of two or more 
jets ("multijets"), with one of the jets misidentified as 
an isolated lepton, and accompanied by large I^t result- 
ing from mismeasurement of jet energies. Significantly 
smaller contributions to the selected sample arise from 
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Z+jets, WW, WZ, ZZ, and single top quark production. 
Together, these five smaller backgrounds are expected to 
contribute from 1% to 7% of the selected sample, depend- 
ing on the number of 6-tagged jets, and are referred to 
below as "other" backgrounds. 

Normalization of the backgrounds begins with the de- 
termination of the number of multijet events in the se- 
lected sample. The multijet background is determined 
using control data samples and probabilities for jets 
to mimic isolated lepton signatures, also derived from 
data 0- Subtracting this background also provides the 
fraction of events with a truly isolated high-p-r lepton (i.e. 

and all backgrounds, except multijets). The contribu- 
tions from single top quark, Z-|-jets, and diboson produc- 
tion are determined from Monte Carlo simulation (MC). 
The remainder corresponds either to ti or VF-|-jet produc- 
tion. The signal and background processes are generated 
using ALPGEN .12] with rut = 175 GeV. pythia is 
used for fragmentation and decay. B hadron decays are 
modeled via evtgen [iJL A full detector simulation is 
performed using GEANT [l5j. 

In an analysis based on the SM, with i? w 1, the tt 
event tagging probabilities are computed assuming that 
each of the signal events contains two 6-jets In the 
present analysis, the top quark can also decay into a 
light quark (d or s) and a W boson. The ratio R de- 
termines the fraction of ti events with 0, 1, and 2 6- 
jets and therefore how ti events are distributed among 
the 0, 1, and 2-tag samples. In order to derive the ti 
event tagging probability as a function of R, we deter- 
mine the tagging probability for the three following sce- 
narios (i) ti W^b W~b (to be referred to as tt bb), 
(ii) ti W^b W^qi or its charge conjugate (referred 
to as tt bqi), and (iii) ti W^qi W~qi (referred to 
as tt qiqi), where qi denotes either a d or s quark. 
The probabilities Pntag to observe ntag = 0, 1, or > 2 
^-tagged jets are computed separately for the three types 
of ti events, using the probabilities for each type of jet 
(&, c, or hght-quark jet) to be 6-tagged. The proba- 
bilities Pntag ill the three scenarios are then combined 
to obtain the ti tagging probability as a function of _R, 

Pn,^M = R^Pn..,{tt ^ bb) + 2R(l - R)Pn,Jtt ^ 

bqi) + {I - R)'^Pnt^g{tt — > qiqi), where the subscript ntag 
runs over 0, 1, and > 2 tags. Table 13 compares the ob- 
served number of events in the 0, 1, and 2-tag samples 
with the sum of the predicted backgrounds and the fitted 
number of ti events. 

The fraction of ti events in the £ + > A jets (^ -I- 3 jets) 
0-tag sample changes from 10% (2%) for i? = 1 to 22% 
(4%) for R = 0. The size of this contribution is of the or- 
der of the Poisson uncertainty on the number of events in 
the 0-tag sample. Therefore the number of observed 0-tag 
events is a poor constraint on R and Nf^. We achieve a 
tighter constraint on the number of ti events in the 0-tag 
sample by constructing a discriminant function V for 0- 
tag events in the £ + >4 jets sample, that combines kine- 



matical event properties to discriminate between ti signal 
and M^-t-jets background. The signal to background ra- 
tio in the ^ -1- 3 jets, 0-tag sample is five times smaller 
than in the corresponding > 4 jets sample. Therefore we 
do not consider such a discriminant for £ + 3 jets, 0-tag 
events. We select four variables that provide good dis- 
crimination between signal and background and that are 
well modeled by the MC. The discriminant function is 
built from: (i) the event sphericity S, constructed from 
the four-momenta of the jets, (ii) the event centrality 
C, defined as the ratio of the scalar sum of the pr of 
the jets to the scalar sum of the energies of the jets, {iii) 
K^^^^ = MlJ^'^p^''' / , where A7^5^'" is the minimum 
separation in — </> space between pairs of jets, p^'" is 
the pt of the lower-pT jet of that pair, and is the 
scalar sum of the lepton transverse momentum and -^t, 
and {iv) Hrp2 — Ht2/Hz, where Ht2 is the scalar sum 
of the Et for all jets excluding the leading jet and 
is the scalar sum of the absolute value of the momenta 
of all the jets, the lepton and the neutrino along the z- 
direction |l7| . Sphericity and centrality characterize the 
event shape and are described in Ref. 0|. In order to 
reduce the dependence on modeling of soft radiation and 
the underlying event, only the four highest-p^ jets are 
used to determine these variables. 

The discriminant function is constructed using the 
method described in Ref. Neglecting correlations 

among the input variables xi,a;2,..., the discriminant 
function can be approximated by the expression: 

n,s,(a;,)A(a;,) + l ' 

where Si{xi) and bi{xi) are the normalized distributions 
of variable Xi for signal and background, respectively. As 
constructed, the discriminant peaks near zero for back- 
ground, and near one for signal. The shapes of the dis- 
criminant for ti and iy-|-jets events are derived from MC. 

The shape of the discriminant for the multijet back- 
ground is obtained from a control data sample, selected 
by requiring that the lepton candidates fail the isolation 
criteria. The other backgrounds (Z-|-jets, diboson, and 
single top quark) have discriminant distributions close to 
those of the VF-|-jet events, and contribute to 1% of the 
0-tag sample. In the final fit, we assume that these pro- 
cesses have the same discriminants as the MK-l-jets events. 
The background normalization in the £+ > A jets, 0-tag 
sample is extracted from the discriminant fit rather than 
from MC. To verify that the kinematic variables used 
in the discriminant are well modeled by the simulation 
we compare data and MC distributions in two control 
samples. To avoid biasing the measurement with respect 
to R, we choose control samples where 6-tagging is not 
applied, and to avoid bias with respect to iVj^ we select 
events with little ti content: £ + 2 jets and ^ -I- 3 jets. In 
£ + 2 jets events, the fraction of ti events is negligible, 
whereas it makes up about 5% of the £ + 3 jets events. 
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In order to measure R and N^t, we perform a binned 
maximum likelihood fit. The data are binned in thirty 
bins: (i) twenty bins of the discriminant 2? in the e + 
>4 jets and ^ + >4 jets, 0-tag samples, (ii) two bins for 
the two 0-tag samples in e + 3 jets and /i + 3 jets, (m) 
four bins for the four 1-tag samples (electron or muon 
and 3 or 4 jets), and (w) four bins for the four 2-tag 
samples (electron or muon and 3 or 4 jets). In each bin, 
we predict the number of events that corresponds to the 
sum of the expected background and signal. The sig- 
nal contribution is a function of R and N^t. To predict 
the number of events in each bin of the discriminant D, 
we use its expected distribution for M^+jets background 
and tt signal. As described earlier, the normalization of 
the multijet background is estimated by counting events 
in orthogonal control samples. Statistical fluctuations in 
the number of events in the control samples are taken 
into account. We incorporate systematic uncertainties 
into the likelihood by using nuisance parameters jl^ . AH 
preselection efficiencies, tagging probabilities, and shapes 
of the discriminant T) are functions of the nuisance pa- 
rameters. The likelihood contains one Gaussian term for 
each nuisance parameter. The value of R that maximizes 
the total likelihood is i? = 1.031° (stat-f syst), m good 
agreement with the SM expectation. A summary of sta- 
tistical and systematic uncertainties is given in Table ITll 
The fit also yields the total number of tt events in the 0, 
1, and 2-tag samples, Nfj = 163^27 (stat). The result of 
the two-dimensional fit is shown in the (i?, Nf^) plane in 
Fig. n^a), with the 68% and 95% contours of statistical 
confidence. In Fig. ^b) and Fig. Q^c), we compare the 
observed number of events to the sum of the predicted 
backgrounds and the fitted tt contribution, in the 0, 1 
and 2-tag samples for events with 3 jets and > 4 jets. In 
Fig. n^d), we compare the observed distribution of the 
discriminant V with the corresponding distribution for 
the sum of the predicted backgrounds and the fitted tt 
contribution. 

We extract lower limits on R and the CKM matrix 
element |Vffc| assuming \Vtb\ = VR- Using a Bayesian 
approach with the prior 7r(i?) = 1 for < i? < 1 and 
7r(i?) = otherwise, we obtain R > 0.78 at the 68% 
C.L. and R > 0.61 at the 95% C.L. For the CKM matrix 
element \Vtb\, we obtain |Vtb| > 0.88 at 68% C.L., and 
\Vtb\ > 0.78 at the 95% C.L. 

In summary, we performed the most accurate measure- 
ment of R to date, R = 1.03to'i7 (stat -|- syst), in good 
agreement with the SM. 
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TABLE I: Observed number of events, predicted backgrounds 
and fitted Ntt. 
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TABLE II: Summary of statistical and systematic uncertain- 
ties on ii. 

Uncertainties on R 

Statistical +0.17 -0.15 

6-tagging efficiency +0.06 —0.05 

Background modeling +0.05 —0.04 

Jet identification and energy calibration +0.04 —0.03 
Multijet background +0.02 

Total error +0.19 -0.17 
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